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Abstract 

The ràpid progress of computer technology has been accompanied by 
a corresponding evolution of software development, from hardwired com- 
ponents and binary machine code to high level programming languages, 
which allowed to master the increasing hardware complexity and fully 
exploit its potential. 

This paper investigates, how classical concepts like hardware abstrac- 
tion, hierarchical programs, data types, memory management, flow of 
control and structured programming can be used in quantum computing. 
The experimental language QCL will be introduced as an example, how 
elements like irreversible functions, local variables and conditional branch- 
ing, which have no direct quantum counterparts, can be implemented, and 
how non-classical features like the reversibility of unitary transformation 
or the non-observability of quantum states can be accounted for within 
the framework of a procedural programming language. 

1 Quantum Programming 

1.1 Quantum Programming Languages 

From a software engineering point of view, we can regard the formalism of 
Hilbert-space àlgebra as a specification language, as the mathematical descrip- 
tion of a quantum algorithm is inherently declarative and provides no means to 
derive a unique decomposition into elementary operations for a given quantum 
hardware. 

Low level formalisms like quantum circuits 0], on the other hand, are usually 
restricted to specific tasks, such as the description of unitary transformations, 
and thus lack the generality to express all aspects of non-classical algorithms. 

The purpose of programming languages is therefore twofold, as they allow 
to express the semàntics of the computation in an abstract manner, as well as 
the automated generation of a sequence of elementary operations to control the 
computing device. 
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1.2 Quantum Algorithms 

In its simplcst form, a quantum algorithm merely consists of a unitary transfor- 
mation and a subsequent measurcmcnt of the resulting state. For more "tradi- 
tional" computational tasks, however, as e.g. searching or mathematical calcu- 
lations, efficient quantum implementations often have the form of probabilistic 
algorithms. Figure ^ shows the bàsic outlinc of a probabilistic non-classical 
algorithm with a simple evaluation loop. 




[ STOP j 

Figure 1: a simple non-classical algorithm 



More complex quantum algorithms, as e.g. Shor's algorithm for quantum 
factoring |SJ |7|, can also include classical random numbers, partial measure- 
ments, nested evaluation loops and múltiple termination conditions; thus the 
actual quantum operations are embedded into a classical flow-control frame- 
work. 

In the discussion of non-classical algorithms, both aspects, the classical con- 
trol structure and the actual quantum operations, are usually treated separately 
and often, only the latter is formally described In order to provide a con- 
sistent formalism, a quantum programming language will have to resolve this 
antagonism by generalizing existing classical programming concepts to the field 
of quantum computing. 



1.3 Structured Quantum Programming 

In traditional computing science, programming languages can be categorized as 
either logical (e.g. Prolog), functional (e.g. LISP) or procedural (e.g. Fortran, 
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Pascal), thc latter bcing the most widely used, for the description of algorithms, 
as well as the actual implementation of real world programs [I]. 
Procedural programming languages can be characterized by 

• explícit flow of control 

• hierarchical program structure 

• tight mapping between code and computation 

which seems to fit most people's way of reasoning about computational tasks. A 
procedural language is called structured, if flow control is restricted to selection- 
and loop-statements with well defined entry- and exit-points (e.g. Modula, 
Pascal without goto-statement) 

Structured quantum programming is about cxtending these concepts into 
the field of quantum computing while preserving their classical semàntics. The 
following table gives an overview of quantum language elements along with their 
classical counterparts. 



classical concept 


quantum analog 


classical machine model 


hybrid quantum architecture 


variables 


quantum registers 


subroutines 


unitary operators 


argument and return types 


quantum data types 


local variables 


scratch registers 


dynamic memory 


scratch space management 


boolean expressions 


quantum conditions 


conditional execution 


conditional operators 


selection 


quantum if-statement 


conditional loops 


quantum forking 


none 


inverse execution of operators 


none 


quantum measurement 



1.4 Hybrid Architecture 

Structured quantum programming uses a classical universal language to definc 
the actual sequence of elementary instructions for a quantum computer, so a 
program is not intended to run on a quantum computer itself, but on a (prob- 
abilistic) classical computer, which in turn controls a quantum computer and 
processes the results of measurements. In the terms of classical computer sci- 
ence, this architecture can be described as a universal computer with a quantum 
oracle (ngure0|. 

From the perspective of the user, a quantum program behaves exactly like 
any other classical program, in the sense that it takes classical input, such as 
startup parameters or interactive data, and produces classical output. 

The state of the controlling computer (i.e. program counter, variable vàlues, 
but also the mapping of quantum registers) is referred to as program state. The 
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í QCL program j 



quantum operations 
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binary 
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quantum 
machine state 



Figure 2: hybrid quantum architecture 



quantum computer itself does not require any control lògic, its computational 
state can therefore be fully described by the common quantum state l'S) of its 
qubits (machine state). 



2 The Programming Language QCL 

QCL (an acronym for "quantum computation language" ) is an experimental 
structured quantum programming language ^O]- A QCL interpreter, written in 
Cl·l — h, including a numerical simulation library (libqc) to emulate the quantum 
backend is available from 

http : //tph . tuwien . ac . at/ ~oe mer/ qcl . html 

as free software under the terms of the GPL. 



2.1 Quantum Storage and Registers 

The smallest unit of quantum storage in QCL is the qubit 

Definition 1 (Qubit) A qubit or quantum bit is a quantum system whose state 
\tp) G B can be fully described by a superposition of two orthonormal eigenstates 
labeled |0) and 1), i.e. B = C 2 . 

|V)=a|Q)+0|l) with |a| 2 + |/3| 2 = l (1) 

QCL treats qubits as quantum registers of length 1. 

qcl> qureg a[l] ; qureg b[l]; // allocate 2 qubits 

qcl> Rot(-pi/3,a) ; // rotate lst qubit 

qcl> H(b) ; // Hadamard Transf ormation 

Definition 2 (Machine State) The machine state \^f) e H of an n-qubit 
quantum computer is the state of a composite system of n identical qubits, i.e. 

n = B® n = c 2 " . 
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QCL — if used together with a numerical simulator — provides debugging 
functions which allow the inspcction of thc othcrwise unobservable machine 
state. In the example below |*) = |00) <g> (H |0}) ® R x (ir/3) |0>: 

qcl> dump; // show product state as generated above 

: STATE: 2/4 qubits allocated, 2/4 qubits free 

0.612372 |0000> + 0.612372 |0010> + 0.353553 |0001> + 0.353553 |0011> 

Definition 3 (Quantum Register) An m qubit quantum register s is a se- 
quence of mutually different qubit positions (sq,si . . . s m _i) of some machine 
state |*) <E C 2 ™ with n>m. 

Using an arbitrary permutation tt over n elements with 7T, = Sj for i < m, a 
unitary reordering operator Tl s is defined as (see also \2.2^\) 

II S \d 0l di . . .d n -i) = Idffojdjn . ■■d 1Tn _ l ) (2) 

Quantum registers are the fundamental quantum data-type in QCL. As they 
contain the mapping between the symbolic quantum variables and the actual 
qubits in the quantum computer, they are the primary interface between the 
classical frontend and the quantum backend of the hardware architecture, as 
the machine state can only be accessed via registers. 

Quantum registers are dynamically allocated and can be local. Temporary 
registers can be created by using the qubit- (q[n]), subregister- (q[n:m]) and 
concatenation-operators (q&p). 



2.2 Operators 

2.2.1 Register Operators 

Definition 4 (Register Operator) The register operator U(s) for an m-qubit 
unitary operator U : C 2 — > C 2 and an m-qubit quantum register s on an n- 
qubit quantum computer is the n-qubit operator 

U(8)=TÚ s (U<g)I(n-m))TL a (3) 

with an reordering operator H s and the k-qubit identity operator I(k). 
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Figure 3: the register operator U(s) 
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All operators in QCL are register operators and can also have an arbitrary 
number of classical parameters. The length of the operand-registers is passed 
as an additional implicit parameter. 

operator dft(qureg q) { // Quantum Fourier Transform 
const n=#q; // set n to length of input 

int i ; int j ; // declare loop counters 

for i=l to n { 

for j=l to i-1 { // apply conditional phase gates (see 2.3.2) 
if q[n-i] and q[n-j] { Phase(pi/2~ (i-j ) ) ; } 

} 

H(q[n-i]); // qubit rotation 

} 

flip(q); // swap bit order of the output 

} 

QCL operators are unitary and have mathcmatical semàntics i.e. their effect 
must be reproduceable and may only depend on the specified parameters. This 
esp. excludes 

• dependències on the program state (e.g. global variables) 

• side effects on the program state 

• user input and calls to randomO 

• non-unitary quantum operations (i.e. calls to measure and reset) 

For any QCL operator, the adjoint operator is dctermined on the fly if the 
call is prefixed with the inversion-flag (I). 1 

qcl> qureg q[2] ; // allocate a 2-qubit register 

qcl> dft(q); // discrete Fourier transform 

[2/32] 0.5 |00> + 0.5 |01> + 0.5 |10> + 0.5 |11> 
qcl> !dft(q); // inverse transform 

[2/32] 1 |00> 

2.2.2 Quantum Data Types 

Classical programming languages often allow to impose access restrictions on 
variables and subroutine parameters. This is equally done to prevent subsequent 
programming errors, as to provide information to the compilcr to allow for more 
cfficicnt optimizations. 

QCL extends this concept to quantum registers by introducing quantum 
data types to limit the ways how operators may effect the machine state. 





type 


restriction 




qureg 


none 




quconst 


invariant to all suboperators 




quvoid 


has to be empty when the uninverted operator is called 




quscratch 


has to be empty before and after the call 



1 Internally, this is achieved by recursively caching all suboperator calls and executing them 
in reverse order with swapped inversion-flags. 
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Definition 5 (Invariance of Registers) A quantum register c is invariant 
to a register operator U(s,c) iffU\i) s \j) c — (Uj \i) s ) \j) c with unitary Uj. 

Definition 6 (Empty Registers) A register e is empty iff i») = |0) e |V') 
2.2.3 Quantum Functions 

Onc important aspect of quantum computing is, that — due to the linearity of 
unitary transformations — an operator applied to a superposition state Y$>) is 
simultaneously applied to all basis vectors that constitute \^>) (quantum paral- 
lelism) since 

^ £ = £ c^i*» ( 4 ) 

i i 

In many cases U implements a reversible boolean, or, equivalently, a bijective 
integer function, by treating the basis vectors merely as bitstrings or binary 
numbers. 

Definition 7 A n-qubit quantum function is a unitary operator of the form 
U : \i) — > \ttí) with some permutation tt over Z 2 «. 

Quantum functions are implemented by the QCL subroutine-type quf unct. 

qufunct inc(qureg x) { // increment register 

int i ; 

for i = #x-l to 1 step -1 { 

CNot (x [i] ,x [0 : i-1] ) ; // apply controlled-not from 
} // MSB to LSB 

Not(x[0]); 

} 

To enforce the above restrictions, the 4 QCL subroutines types form a calling 
hierarchy, i.e. routines may only invoke subroutines of the same or a lower level. 



subroutine 


s 


* 


invertible 


description 


procedure 


all 


all 


no 


classical control structure 


operator 


none 


unitary 


yes 


general unitary operator 


qufunct 


none 


permutation 


yes 


quantum function 


functions 


none 


none 


no 


mathematical functions 



The columns S and tp denote the allowed side-efïects on the classical program 
state and the quantum machine state. 



2.2.4 Irreversible Functions 

One obvious problem in QC is its restriction to reversible computations. Con- 
sider a simple operation like computing the parity of a bitstring 

parity' : \i) -> |6(i)mod2) with b(n) = n mod 2 + b( [n/2\ ), 6(0) = (5) 

Clearly, this operation is non-reversible since parity' |1) = parity' |2), so 
parity' is not an unitary operator. However, if we use an additional target 
register, then we can definc an operator parity which matches the condition 
parity |í,0) = \i, b(i) mod 2). 



7 



qufunct parity (quconst x, quvoid y) { 
int i ; 

for i = to #x-l { 

CNot(y,x[i] ) ; // flip parity for each set bit 

> 

} 

In QCL, an operator F : |a;) x |0) y — * \x) x \f(x)) y is declared as qufunct with 
at least one invariant (quconst) argument register x and one empty (quvoid) 
target register y as parameter. The result for -F'|a;) x |y 7^ 0) y is unspecified to 
allow for different ways to aceumulate the result. So F 1 : \x, y) — > \x,y® f{x)) 
and F" : \x,y) — » \x + y mod 2 k ) are merely considered to be different imple- 
mentations of the same quantum function. 

2.2.5 Scratch Space Management 

While quantum functions can be used to work around the reversible nature 
of QC, the necessity to keep a copy of the argument is a problem, as longer 
computations will leave registers filled with intermediate results. 

Let F be a quantum function with the argument register x (quconst), the 
target register y (quvoid) and the scratch register s (quscratcli) 

^(x ) y !S ):|i) x |0) y |0) s ^|i)J/(i)) y |i(i)) s (6) 

F fills the register s with the temporary junk bits To reclaim s, QCL 

transparently allocates an auxiliary register t and translates F into an operator 
F' which is defined as 

F'(x,y,s,t) = J Pt( X) t,s)fanout(t,y)F(x,t,s) (7) 

The fanout operator is a quantum function defined as 

fanout : |i)|0) \i)\i) (8) 

The application of F' restores the scratch register s and the auxiliary register t 
to |0) while preserving the function value in the target register y: 

|i, 0,0,0) -► |i,0 > j(i)>/(0) |t\/(i)JW,/(<)> -> !*./(*). 0,0) (9) 

2.3 Quantum Conditions 
2.3.1 Conditional Operators 

Classical programs allow the conditional execution of code in dependence on the 
content of a boolean variable (conditional branching) . 

A unitary operator, on the other hand, is static and has no infernal flow- 
control. Nevertheless, we can conditionally apply an n-qubit operator U to a 
quantum register by using an enable qubit and define an n + 1 qubit operator 

"-(V u) m 
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So U is only applied to base- vectors where the enable bit is set. This can be 
easily extended to enable-registers of arbitrary length. 



Definition 8 A conditional operator f7[[ e ]] with the enable register e is a unitary 
operator of the form 



U M] :\i,e) = \i)\e) e ^í. 




A conditional version of the increment operator from l2.2"3l can be explicitly 
implemented as 

qufunct cinc(qureg x,quconst e) { 



QCL can automatically derive í/[[ e ]] for an operator U if its declaration is 
prefixed with cond. 

cond qufunct inc(qureg x.quconst e) { ... } 

The enable register can be set by a quantum if-statement and QCL will trans- 
parently transform the defined operator into its conditional version when nec- 
essary. 

2.3.2 Quantum if-statement 

Just like the concept of quantum functions allows the computation of irreversible 
boolean functions with unitary operators, conditional operators allow condi- 
tional branching depending on unobservable qubits. 

Given the above defmitions, the statement if e { inc (x) ; } is equivalent 
to the explícit call of cinc (x,e) and if e { inc(x); } else { !inc(x); } 
is equivalent to the sequence 

cinc(x,e); // conditional increment 

Not(e); // invert enable qubit 

!cinc(x,e); // conditional decrement 

Not(e); // restore enable qubit 

Quantum if-statements can be nested. Since operators within the if- and else- 
branches are transformed into their conditional versions, they must be declared 
cond and must not operate on any qubits used in the condition. 

2.3.3 Complex Conditions 

Conditions in quantum if-statements are not restricted to single qubits, but 
can contain any boolean expression and also allow the mixing of classical and 
quantum bits. 



int i ; 

for i = #x-l to 1 step -1 { CNot (x [i] ,x [0 : i-1] & e) ; } 
CNot(x[0] ,e) ; 



} 
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qcl> qureg q[4] ; qureg b[l]; qureg a[l] ; 
qcl> H(a & b) ; // prepare test-state 

[6/32] 0.5 |000000> +0.5 |010000> +0.5 | 100000> +0.5 | 110000> 
qcl> if a and b { inc(q); } 

[6/32] 0.5 |000000> + 0.5 |010000> + 0.5 | 100000> + 0.5 |110001> 
qcl> if a or b { inc(q); } 

[6/32] 0.5 |000000> + 0.5 |010001> + 0.5 | 100001> + 0.5 |110010> 
qcl> if not (a or b) { inc(q); } 

[6/32] 0.5 |000001> +0.5 |010001> +0.5 | 100001> +0.5 | 110010> 

QCL produces a sequence of CNot-gates to evaluate a quantum condition. 2 
If necessary, scratch qubits are transparently allocated and uncomputed again. 



2.3.4 Forking if-statement 

If the body of a quantum if-statement contains statements which change the pro- 
gram state (e.g. assignments to local variables), then subsequent operator calls 
may differ, depending on whether the if- or the else-branch has been executed. 

In that case, QCL follows all possible classical paths throughout the op- 
erator definition (forking), accumulates the conditions of all visited quantum 
if-statements and serializes the generated sequence of operators. 

cond qufunct demux(quconst s, qureg q) { 
int i ; 
int n = 0; 

for i=0 to #s-l { // accumulate content of 

if s[i] { n=n+2~i; } // selection register in a 

} // classical variable 

Not(q[n]); // flip selected output qubit 

} 

Figure^shows the quantum circuit [Sj generated by demux(s,q) in the case 
of a 2-qubit selection register s. 
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Figure 4: a quantum demultiplexer 



Forking if-statements may only appear within an operator definition to as- 
sure that the different execution threads can be joined again. 

2 Internally, this is achieved by transforming a quantum condition into its exclusive dis- 
junctive normal form 
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3 Conclusions 



Throughout the history of classical computing, hardware development has al- 
ways been accompanied by a corresponding improvcmcnt of programming method- 
ology. The formalism of quantum circuits — the moral equivalent to hand writ- 
ten machine code — seems inadequate for quantum computers with more than 
a couple of qubits, so more abstract methods will cvcntually be required. 

We have demonstrated how wcll cstablished concepts of classical program- 
ming languages like subroutines, local variables or conditional branching can 
be ported to the ficld of quantum computing. Besides providing a new level 
of abstraction, we also hope that a quantum programming language which se- 
mantically integrates those concepts will allow for a better and more intuitive 
understanding of non-classical algorithms. 
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